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Abstract 

In recent twenty years, loop quantum gravity, a background independent ap- 
proach to unify general relativity and quantum mechanics, has been widely inves- 
tigated. We consider the quantum dynamics of a real massless scalar field coupled 
to gravity in this framework. A Hamiltonian operator for the scalar field can be 
well defined in the coupled diffeomorphism invariant Hilbert space, which is both 
self-adjoint and positive. On the other hand, the Hamiltonian constraint operator 
for the scalar field coupled to gravity can be well defined in the coupled kinemat- 
ical Hilbert space. There are 1 -parameter ambiguities due to scalar field in the 
construction of both operators. The results heighten our confidence that there is 
no divergence within this background independent and diffeomorphism invariant 
quantization approach of matter coupled to gravity. Moreover, to avoid possible 
quantum anomaly, the master constraint programme can be carried out in this cou- 
pled system by employing a self-adjoint master constraint operator on the diffeo- 
morphism invariant Hilbert space. 

Keywords: loop quantum gravity, quantum dynamics, scalar field, background in- 
dependence. 
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1 Introduction 

The research on quantum gravity theory is rather active. Many quantization pro- 
grammes for gravity are being carried out (for a summary see e.g. Q]). In these differ- 

*Email address: mhanl@lsu.edu 
^Email address: mayg@bnu.edu.cn 



1 



ent kinds of approaches, the idea of loop quantum gravity is motivated by researchers 
in the community of general relativity. It follows closely the thoughts of general rela- 
tivity, and hence it is a quantum theory born with background independency For 
recent review in this field, we refer to 1 3 1 1 4 1 1 5 1 . 

To apply the background independent quantization technique, one first casts gen- 
eral relativity into the Hamiltonian formalism of a diffeomorphism invariant Yang- 
Mills gauge field theory with a compact internal gauge group [6||7|. The kinematical 
Hilbert space of the quantum theory is then constructed rigorously. One can even 
solve the Gaussian and diffeomorphism constraints to arrive at a diffeomorphism in- 
variant Hilbert space |8|. Certain geometrical operators are shown to have discrete 
spectra in the kinematical Hilbert space |9||10||11||12||13|. However, some impor- 
tant elements in this approach are not yet understood. Despite the systematic efforts 
in constructing the Hamiltonian constraint operatorl 14|| 15 1 and the master constraint 
operatorl 16||17||18|, the dynamics of the quantum theory has not been fully under- 
stood, especially if one wants to include all known matter fields. The primary goal 
of this paper is to apply the loop quantization technique to a scalar field and check 
whether the quantum dynamics can be well defined. Since we will use the developed 
polymer-like kinematical description of the scalar field| 19 1 1 20 1, it can be considered as 
a development of the construction for quantum Higgs fields in Refs. |21 1 and |22|. 

In section 2, the Hamiltonian formalism of a massless real scalar field coupled to 
gravity is obtained in generalized Palatini formulation. For readers' convenience, the 
loop quantum kinematical setting of a real scalar field coupled to gravity is also intro- 
duced. We then show in section 3 that an operator corresponding to the Hamiltonian 
of the scalar field can be well defined on the coupled diffeomorphism invariant Hilbert 
space. It is even positive and self-adjoint. Thus quantum gravity acts exactly as a nat- 
ural regulator for the quantum scalar field in the polymer representation. In section 4, 
to study the whole dynamical system of the scalar field coupled to gravity, a Hamilto- 
nian constraint operator is defined in the coupled kinematical Hilbert space. Moreover, 
the contribution of the scalar field to the Hamiltonian constraint can be promoted to a 
positive self-adjoint operator. Similar to the gravitational Hamiltonian constraint, there 
is an one-parameter ambiguity in defining both the Hamiltonian operator and the con- 
straint operator due to the scalar field. To avoid possible quantum anomaly and find 
physical Hilbert space, the programme of master constraint for the coupled system is 
discussed in section 5. A self-adjoint master operator is obtained in the diffeomorphism 
invariant Hilbert space, which assures the feasibility of the programme. 

2 Polymer-like representation for scalar field coupled 
to gravity 

Consider first the classical dynamics of a real massless scale field (f> coupled to gravity 
on a 4-manifold M. The coupled generalized Palatini action reads 

S[4,oj^,<p] = 5 p [4,^ /y ] +S KG [e fS K ,<f>l 
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where 



Z JM 

here e^, and oj^ j are respectively the tetrad and Lorentz connection on M, the real 
number y, k and % are respectively the Barbero-Immirzi parameter, the gravitational 
constant and the coupling constant. From now on we use a, /?,••• for 4-dimensional 
spacetime indices and /, J, ■ ■ ■ for internal Lorentz indices, a, b, ■ ■ ■ for 3-dimensional 
spacial indices and i, j, ■ ■ ■ for internal SU(2) indices. After 3+1 decomposition and 
Legendre transformation, similar to the case in Palatini formalism|23|, we obtain the 
total Hamiltonian of the coupling system on the 3-manifold X as: 

<H M = f(A i G i +N a <V a +NC), (1) 

where A', N" and N are Lagrange multipliers, and the Gaussian, diffeomorphism and 
Hamiltonian constraints are expressed respectively as: 

d = D a I» := d a P? + e ij k A i a P a k , (2) 
= PtF^-A^Gt+W^, (3) 

C = K 2— P c lP h j [e li k F k ab - 2(1 + y 2 )K[ a Ki ] 



_1 f^^^p b{dam<j) + 1 ^ 
Vldet^l 2 2 «m 



(4) 



here the conjugate pair for gravity consists of the SU(2) connection A' a and the densi- 
tized triad P b ., F' ab is the curvature of A' a , and n denotes the momentum conjugate to (p. 
Thus one has the elementary Poisson brackets 

{AUxXP'/y)} = 5 a b 5)5{x-y), 
{(p(x),n(y)} = 5(x-y). 

Note that the second term of the Hamiltonian constraint (0} is just the Hamiltonian of 
the real scalar field. 

Now we introduce the background independent quantization of a real scalar field 
coupled to gravity, following the polymer representation of the scalar field [19|[20|. 
The classical configuration space, 1A, consists of all real-valued smooth functions <p on 
S. Given a set of finite number of points X = {x\, x^} in S, denote Cylx the vector 
space generated by finite linear combinations of the following functions of </>: 
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where A = (A\, A2, ■ ■ •, An) are arbitrary real numbers. It is obvious that Cylx has the 
structure of a *-algebra. The space Cyl of all cylindrical functions on H is defined by 

Cyl := UxCylx. (5) 

Completing Cyl with respect to the sup norm, one obtains a unital Abelian C*-algebra 
Cyl. Thus one can use the GNS structure to construct its cyclic representations. A 
preferred positive linear functional lqq on Cyl is defined by 

/TT , f 1 iU» = 0Vj 

m{n x , A ) = I Q oth ^ mise) 

which defines a diffeomorphism-invariant faithful Borel measure fj. on 1/ as 



f 



otherwise. (6) 



Thus one obtains the Hilbert space, < Hf^ = L 2 (U, d/-i), of square integrable functions 
on a compact topological space 1A with respect to fi, where Cyl acts by multiplication. 
The quantum configuration space U is called the Gel'fand spectrum of Cyl. More 
concretely, for a single point set Xq = \xq\, Cylx is the space of all almost periodic 
functions on a real line R. The Gel'fand spectrum of the corresponding C*-algebra 
Cyl Xg is the Bohr completion R Ao of R 1 19 1, which is a compact topological space such 
that Cyl Xg is the C*-algebra of all continuous functions on R Co . Since R is densely 
embedded in R Ao , R A|) can be regarded as a completion of R. 

Given a pair (xq, Aq), there is an elementary configuration for the scalar field, the 
so-called point holonomy, 



U(x ,A ) := exp[iA <p(x )]. 
It corresponds to a configuration operator U(xo, Aq), which acts on any cylindrical func- 

JKC 

'kin 



tion ^(0) 6 ^1° by 



U(x , A M<f>) = U(x , Ao)H<P)- (7) 

All these operators are unitary. But since the family of operators U(xq, A) fails to be 
weakly continuous in A, there is no operator 0(x) on < H^S. The momentum functional 
smeared on a 3-dimensional region R c E is expressed by 

n(R) := I d 3 xn(x). 
Jr 

The Poisson bracket between the momentum functional and a point holonomy can be 
easily calculated to be 

\7t(R), U(x, A)} = -iA X R(x)U(x, A), 

where xr(x) is the characteristic function for the region R. So the momentum operator 
is defined by the action on scalar network functions H c= (x,A) as 

7t{R) II C <#) := m{n(R),Il c ((f>)} = fifj] ^ fi (x,-)]II c (#). 
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It is clear from Eq.© that an orthonormal basis in 1~lE G is given by the so-called scalar 
network functions n c (0), where c denotes (X{c), A) and A = (A\, A%, ■ ■ ■, A^) are non- 
zero real numbers now. So the total kinematical Hilbert space Jitm is the direct product 
of the kinematical Hilbert space for gravity and the kinematical Hilbert space 
for real scalar field, i.e., < Hki n '■= ® ^kin ■ ^ et n. s (A) be the spin network basis 
in <Hgs labeled by s (TT]lE3l- Then the state n s , c = Tl s (A) ® n e (0) e Cyl y(s) (lA) ® 
Cylx( C ) = Cyly( StC )(^l) ® Cyly( StC )(1() is a gravity-scalar cylindrical function on graph 
y(s, c) = y(s) UX(c). Note that generally X{c) may not coincide with the vertices of the 
graph y(s). It is straightforward to see that all of these functions constitutes a complete 
set of orthonormal basis in < Kj ; ,„ as 

< n sV ,(A,0)|n s>c (A,^) > ktn = 6 S > S 6 C > C . 

Note that none of < 74,„, 'H^ and is a separable Hilbert space. 

Now we can consider the quantum dynamics and impose the quantum constraints 
on 'H/tin- Firstly, the Gaussian constraint can be solved independently of "TfEf, since 
it only involves gravitational field. It is also expected that the diffeomorphism con- 
straint can be implemented by the group averaging strategy in the similar way as that 
in the case of pure gravity. Given a spatial diffeomorphism transformation tp, a unitary 
transformation is induced by ip in the Hilbert space 'K / t,„, which is expressed as 

Uq,Tls=(y(s) j,N),e=(X(cM) = ~tt<pas=(<p(y(.s)) j,N)#°c=(p(X(c)Vt) ■ 

Then the differomorphism invariant spin-scalar-network functions are defined by group 
averaging as 

n M := —J— Yj Tj ^^' n ^ ( 8 > 

where Diff y is the set of diffeomorphisms leaving the colored graph y invariant, GS y 
denotes the graph symmetry quotient group Diff y /TDiff y where TDiff Y is the set of 
the diffeomorphisms which is trivial on the graph y, and n y is the number of elements in 
GS y . Following the standard strategy in quantization of pure gravity, an inner product 
can be defined on the vector space spanned by the diffeomorphism invariant spin-scalar- 
network functions such that they form an orthonormal basis as: 

< n[ i>c ]|n[ s / jC /] >Diff'= Tl[s,c][Rs' ,c'£{s' ,c']] - <5[.5,c], [.?',£-']■ (9) 

After the completion procedure, we obtain the expected Hilbert space of diffeomor- 
phism invariant states for the scalar field coupled to gravity, which is denoted by Hoiff- 
In the following sections, we would like to discuss the quantum dynamical properties 
of the polymer-like scalar field coupled to gravity. 

3 Diffeormorphism invariant quantum Hamiltonian of 
scalar field 

In the present section, we first consider the quantum scalar field on a fluctuating back- 
ground. A similar idea was considered in Ref.|25|, where a Hamiltonian operator with 
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respect to a U(l) group representation of scalar field is defined on a kinematical Hilbert 
space Him of matter coupled to gravity. Then an effective Hamiltonian operator of the 
scalar field can be constructed as a quadratic form via 

< ^matter, H n iftter( m ) ^matter > kin' 
■= < 4>gra V (m) ® if/matter, ^matter 4>grav(m) ® ^/ matter >kin>, (10) 

where il/ grm >(m) e H?. R presents a semiclassical state of gravity approximating some 
classical spacetime background m where the quantum scalar field lives. Thus the ef- 
fective Hamiltonian operator H e J^ tter {m) of scalar field contains also the information of 
the fluctuating background metric. In the light of this idea, we will construct a Hamil- 
tonian operator Hkg f° r scalar field in the polymer-like representation. It turns out 
that the Hamiltonian operator can be defined in the Hilbert space Hoijf of diffeomor- 
phism invariant states for scalar field coupled to gravity without UV-divergence. So the 
quantum dynamics of the scalar field is obtained in a diffeomorphism invariant way, 
which is expected in the programme of loop quantum gravity. Thus, here an effec- 
tive Hamiltonian operator of the scalar field could be extracted in Hoijf by defining 
< V F[„,](A, <p), Hkg ^[m] (A, <p) >Diff to be its expectation value on diffeomorphism in- 
variant states ^((p) of scalar field, where the diffeomorphism invariant semiclassical 
state *P[„,](A) represents certain fluctuating geometry with spatial diffeomorphism in- 
variance, and the label [to] denotes the classical geometry approximated by *P[ m ](A). 
Moreover, the quadratic properties of the scalar field Hamiltonian will provide power- 
ful functional analytic tools in the quantization procedure, such that the self-adjointness 
of the Hamiltonian operator can be proved by the theorem in functional analysis. 



3.1 Regularization of the Hamiltonian 

In a suitable gauge, the classical Hamiltonian for the massless scalar field on the 3- 
manifold S can be given by 

Hkg - Hkg,<p + Hkg, Kin 



d 3 x[ ' ^ =s l Jp?p°.(d a cf>)d b <l> + = ; *ri (ID 

i 2 yj\detq\ 2a M y|det^| 

We will employ the following identities: 

pf = ^r^^A and 4oo = ^{K(x),vuj, 

where rj abc denotes the Levi-Civita tensor tensity and Vu x is the volume of an arbitrary 
neighborhood U x containing the point x. By using the point-splitting strategy, the 
regulated version of the Hamiltonian dl 11 is obtained as: 

Hkg,* = f d 3 y f d 3 x Xe (x-y)^x 



' P*(x)(d a <p(x))^—P b J (y)d b <l>(y) 



6 




£d 3 y£d 3 x Xe (x-y)6 ij x 



77 flcc (5 fl 0(x))Tr(T,{A c (x), V^ 4 }{A c (x), V^ 4 }) x 
Tf fd (d b< p(ymirj{A f (y), V 3 J?}{A d (y), vffl), 




J^d 3 wTf ef Tr({A d (w), ^}{A e (w), yJV^ w }{A f (w), yjV^}) x 




d\lj ahc Tr({A a (u), yjV^}{A b (u), VV^){A C ( M ), y/v^}) x 



where we denote A a = A' a Tt, Xe(x - y) is the characteristic function of a box containing 
x with scale e such that \\m e ^oXei x ~ y)/ e3 = S(x ~ y)> an d Vu% is the volume of the 
box. In order to quantize the Hamiltonian as a well-defined operator in polymer-like 
representation, we have to express the classical formula of Hkg m terms of elementary 
variables with clear quantum analogs. This can be realized by introducing a triangu- 
lation T(e) of E, where the parameter e describes how fine the triangulation is. The 
quantity regulated on the triangulation is required to have correct limit when e — > 0. 
Given a tetrahedron A € T(e), we use {5 , ,(A)} (= i j 2,3 to denote the three outgoing oriented 
segments in A with a common beginning point v(A) = s(s,'(A)) and use a,y(A) to denote 
the arcs connecting the end points of s,-(A) and sj(A). Then several loops ajj(A) are 
formed by a,y(A) := s,(A) o a^A) o s 7 (A) _1 . Thus we have the identities: 



for nonzero A, where s(s(A)) and t(s(AJ) denote respectively the beginning and end 
points of segment s(A) with scale e associated with a tetrahedron A. Regulated on the 
triangulation, the classical Hamiltonian of scalar field reads 



{ dt A a s a (t),V 3 ' 4 ) = -A(s(A))- l {A( S (A)),Vl 




) + o(e), 



and 




KG 4 ~ 



4«M 

9*V 



A'er(e) Aer(e) 




Tr(r,A(^(A)r 1 {A(. m (A)),y^ 4 }A(s n (A)y l {A( Sn (A)), V 3 ' 4 }) x 
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e kpq \u(v(A'), A)^ [U(t(s k (A')), A) - U(v(A'),A)] x 



A 



Tr(T J A(. p (A'))- I {A( Sp (A')), V%\ }A{s q {A')T x {A{s q {A')), V%\ }), 



KG ' Km 8lQ- M (^y) 6 ^ ^ 



J] e i ™Tr(A(5 ! (A"))- 1 {A(5 i (A")), ^ 



u e ) x 



A"eT(e) 



A(s m (A")r\A(s m {A")), x 
A(s n {A")y l {A{s n {A")), ^~)) x 
X e-'' w Tr(A(. S /A"'))- 1 {A(^(A"')), x 

A"'€7Xe) 

A(^(A"'))- I {A(^(A"')), A /VVT) x 



<A"' ) 



A( i/ (A"'))- 1 {A(^(A"')), ^~ } ))x 

xJv(A) - v{A')) Xe (y(A") - v(A))^(v(A'") - v(A')), (12) 

where the overall numerical factors are got from the scale of a tetrahedron to an oc- 
tahedron. Note that the above regularization method is in complete analogy with that 
used for Higgs field in Ref. 1 22 1. However, our regularization formula of scalar field 
Hamiltonian is explicitly dependent on the parameter A, which will leads to a kind of 
quantization ambiguity of the real scalar field dynamics in polymer-like representation. 

3.2 Quantization of the Hamiltonian 

Since all constituents in the expression (I12> have clear quantum analogs, one can quan- 
tize it as an operator by replacing the constituents by the corresponding operators and 
Poisson brackets by canonical commutators. Then the regulator should be removed by 
e — > with respect to a suitable operator topology. Now we begin to construct the 
Hamiltonian operator. Given a spin-scalar-network function in order to ensure 
that the final operator is diffeomorphism covariant, and cylindrically consistent (up to 
diffeomorphisms), one can make the triangulation T(e) adapted to the graph y(s, c) of 
n j c according to the strategy developed in Ref.[ 14 1 with the following properties. 

• The graph j(s, c) is embedded in T(e) for all e, so that every vertex v of y(s, c) 
coincides with a vertex v(A) in T(e). 

• For every triple of edges (ei, <?2, ei) of y(s, c) such that v = s(ei) = sfe) = sfe), 
there is a tetrahedra A e T(e) such that v = v(A) and s,(A) c e,-, V i = 1, 2, 3. We 
denote such a tetrahedra as A" Si . 

• For each tetrahedra A^ s s one can construct seven additional tetrahedron A^, p p , 



1, 7, by backward analytic extensions of S;(A) so that t/ JljS2jJ3 := U^A* 3 



1 •■ , 2'"3 

is a neighborhood of v. 
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The triangulation must be fine enough so that the neighborhoods U(v) := U S] S2 J3 t/ JljJ2>J3 (v) 
are disjoint for different vertices v and v' of y(s, c). Thus for any open neigh- 
borhood Uy(^ C ) of the graph y(s, c), there exists a triangulation Tie) such that 

^veV(y(s,c))U(v) C I/y( JiC ). 

The distance between a vertex v(A) and the corresponding arcs a,/(A) is described 
by the parameter e. For any two different e and e', the arcs aij(A € ) and fly(A e ) 
with respect to one vertex v(A) are analytically diffeomorphic with each other. 

With the triangulation 7\e), the integral over S is replaced by the Riemanian 
sum: 

/ - / + l ' 
f = Z f + f 

where n(v) is the valence of the vertex v = s(s\) = s(s2) = s(si,), and C\ = iy ^ 
Defficient. One then o 

f = 8 f 

Jft, ,„,,ftO J A?. 



denotes the binomial coefficient. One then observes that 



in the limit e — > 0. 
• The triangulation for the regions 

f/(v) - U W3 (v), 

Uy( s ,c) ~ U v el/(y(s,e)) C/(v), 

E-EW)' (13) 

are arbitrary. These regions do not contribute to the construction of the oper- 
ator, since the commutator terms like [A(s,(A)), Vt/, i(A) ]n 5jC would vanish for all 
tetrahedron A in the regions Jl 31 . 

Note that there are many possible ways in choosing a triangulation. The motivation 
for above choice has been fully discussed in Ref.| 14 1. Introducing a partition f of the 
3-manifold E into cells C, we can smear the "square roots" of c H e KC< p and 1~i e KC Kin in 
one cell C respectively and promote them as regulated operators in 'Hkin with respect 
to the state-dependent triangulation T(e) as 



w e - c - V f \ V _P A _fe,A Pa 

'Ms,c)) v(A)=v V-^W V fi W 



w: 



veV(y(j,c)) v(A)=v 



Z Z -^r^v-^=. 



v£V(y(i,c)) v(A)=v yzyv) V^W 
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where Xc( v ) is the characteristic function of the cell C, and 

Ka : = ^e"""^^(v^(v))- 1 [f/(f(5,(A))^(v))-j7(v^(v))]x 

Tr(r,A( im (A))- 1 [A( im (A)), V 3 J f 4 ]A(s n (A)y l [A(s„(A)), V^ 4 ]), 
: = ^3^v)e* m Tr(i(s,(A))- I [i(s / (A)), ^] x 
A(^(A))- 1 [A(* m (A)), 

^(Air'tiWA)), a/^]). (15) 
Note that the tetrahedron projector p& associated with segments s\, S2 and 53 reads 

PA := P Sl P S2 Ps 3 

R 1 /l 1 /l 1 

= e(^--A Sl --)e(^--A S2 --)6(^--A S3 --) (16) 

where A s . is the Casimir operator associated with the segment s, and 9 is the distribution 
on R which vanishes on (-00, 0] and equals 1 on (0, 00), and the tetrahedron projector 
related to a same vertex constitute the vertex operator E(v) := Dv(A)=v^A- Note also 
that the partition P is not required to coincide with the triangulation Tie). We have 
arranged the operator p&/ ^E(v) in such a way that both operators in J14i and their 
adjoint operators are cylindrically consistent up to diffeomorphisms. Thus there are 
two densely defined operators W? t and W^ in in < Hkm associated with the two consistent 
families of (I14i . We now give several remarks on their properties. 

• Removal of regulator e 

It is not difficult to see that the action of the operator W^ c) , i on a spin-scalar- 
network function n s c is graph-changing. It adds finite number of vertices with 
representation A(v) at f(s,(A)) with distance e from the vertex v. Recall that 
the action of the gravitational Hamiltonian constraint operator on a spin network 
function is also graph-changing. As a result, the family of operators $^ JC ) ^ , also 
fails to be weakly convergent when e — > 0. However, due to the diffeomorphism 
covariant properties of the triangulation, the limit operator can be well defined 
via the uniform Rovelli-Smolin topology, or equivalently, the operator can be 
dually defined on diffeomorphism invariant states. But the dual operator cannot 
leave Irioiff invariant. 

• Quantization ambiguity 

As a main difference of the dynamics in polymer-like representation from that 
in U(l) group representation \ 22\, a continuous label A appears explicitly in the 
expression of dl4l> . Hence there is an one-parameter quantization ambiguity due 
to the real scalar field. Recall that the construction of gravitational Hamiltonian 
constraint operator also has a similar ambiguity due to the choice of the repre- 
sentations j of the edges added by its action. A related quantization ambiguity 
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also appears in the dynamics of loop quantum cosmology |26|. We will come 
back to this point in a future publication 1271 . Since the regulator is removed 
in a diffeomorphism invariant way, the quadratic form, which we are going to 
constructed, will be independent of the initial triangulation T in the sense that it 
depends only on the diffeomorphism class of T, as is the case of the gravitational 
Hamiltonian constraint operator [ 14 1. 

Since our quantum field theory is expected to be diffeomorphism invariant, we would 
like to define the Hamiltonian operator of polymer scalar field in the diffeomorphism 
invariant Hilbert space For this purpose we fix the parameter A to be a non- 

zero constant at every point. Then what we will do is to employ the new quantization 
strategy developed in Refs. 1 16 1 and 1 17 1. We first construct a quadratic form in the 
light of a new inner product defined in Ref. 11171 on the algebraic dual D* of the space 
of cylindrical functions. Then we prove that the quadratic form is closed. Note that, 
although the calculation by employing this inner product is formal, it can led to a well- 
defined expression of the desired quadratic form Ea. (l22> . Since an arbitrary element of 
D* is of the form *P = £ JjC c s , c < n s c | ■ >u„, one can formally define an inner product 
< ■ |- >* on D* via 



where the Cantor aleph N denotes the cardinal of the set [s, c]. Note that we exchange 
the coefficients on which the complex conjugate was taken in Ref.| 17 1, so that the inner 
product < ^mffV^Diff >* reduces to < ^DiffW Biff >aff for an Y ^Diff,^" Diff e 
'Huif s- Completing the quotient with respect to the null vectors by this inner product, 
one gets a Hilbert space Our purpose is to construct a quadratic form associated to 
some positive and symmetric operator in analogy with the classical expression of d!2i . 
So the quadratic form should first be given in a positive and symmetric version. It is 
then natural to define two quadratic forms on a dense subset of 'Hoiff C "H* as: 




(17) 





(18) 



where the dual limit operator W' c of either family of W^'f or W^f n in J14i is naturally 
defined on diffeomorphism invariant states as 



w' c v Diff [n s , c ] = iimy Diff [w e ' c n s , c ]. 



(19) 
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To show that the quadratic forms are well defined, we write 

%^m ff = 2 < i>s , e (Q < n s , c | • >* => wj x jc) = (W'^DiffW^, 

s,c 
s,c 

Then, by using the inner product < I17I > the quadratic forms in d 1 81 become 

Qm^Diff^mff) 

- 



QKC,Kin^Diff^' Diff ) 

: = lim V 8 4 x ^— V wl„ , r (C)£ (C) 

- 2 84 >< Z xTiTTrt Z ^AQ^Jq. 

(20) 



81a M (,r) 6 ^N([,,c])^ e] 



Note that, since ^Diff is a finite linear combination of the diffeomorphism invariant 
spin-scalar-network basis, taking account of the operational property of W' c there are 
only finite number of terms in the summation D[ s c ] contributing to ( 1201 . Hence we can 
interchange 2[j C ] an d limp->£ 2ce?> m above calculation. Moreover, for a sufficiently 
fine partition such that each cell contains at most one vertex, the sum over cells there- 
fore reduces to finite terms with respect to the vertices of y(s, c). So we can interchange 
Zs,ce[>,c] an d lim*>-»E 2cep to obtain: 

- Mx w^siRi Z feZ-L^)^^ 

/ [s,c] VL ' Jy s.cefoc] CeP 
' [j,c] VL ' Jy s.cefoc] veV(y(j,c)) 

- 84 * ^ Z ^ Z^^f^^K^^tn,,,, 

(21) 
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where the limit f — > S has been taken so that C — > v. Since given 7(5, c) and 7(5', c') 
which are different up to a diffeomorphism transformation, there is always a diffeo- 
morphism tp transforming the graph associated with W^ f) n jc (v e 7(5, c)) to that of 

^', e <) n ^' ( v ' e r(*'.c')) with <p(v) = v', (W' v y Diff )[n s , c£[s , c] ] is constant for differ- 
ent (s, c) G [s, c], i.e., all the N([s, c]) terms in the sum over (s, c) e [s, c] are identical. 
Hence the final expressions of the two quadratic forms can be written as: 



[s,c] veV(y(s,c)) 



QKG,KI«WDiff,y Diff ) 



3 x 



Ml " [j,c] veV(y(j,e)) 



(22) 

Note that both quadratic forms in (I22t have finite results and hence their form domains 
are dense in 'Hay/- Moreover, both of them are obviously positive, and the following 
theorem will demonstrate their closedness. 

Theorem: Both Qkg.ip an d Qkckw o-re densely defined, positive and closed quadratic 
forms on 'Hoiff, which are associated uniquely with two positive self-adjoint operators 
respectively on 'Hmff such that 

QKcA^D if f^' Diff ) = <T / ,„,'/Y A( ,,M'; ) „, >„,,, 

QKCKmCPoiff^Diff) = < X PDiff\'f(KG,Kin\^' Diff >Diff ■ 

Therefore the Hamiltonian operator 

HkG '■- ^KCip + f^KG.Kin (23) 

is positive and also have a unique self-adjoint extension. 

Proof: We follow the strategy developed in Refs. 11171 and HI 81 to prove that both 
Qkg$ an d Q.KG,Kin are closeable and uniquely induce two positive self-adjoint oper- 
ators < Ka:c.(/. and 'Hkg,Kw- One can formally define 'Hkg.q an d ^KG.Km acting on dif- 
feomorphism invariant spin-scalar network functions via: 



"f^KCip ntn.e,] ^ CA , G,0(n[j 2 .c2]'n[si,ci])n[s2,c2]> (24) 

'HKG,Kin II^^cj] := ^ QKG,Kin{^[s2,C2\, n[j 1]Cl ])II[j 2>C2 ] . (25) 



Then we need to show that both of the above operators are densely defined on the 
Hilbert space "Hoi/fi i- e -> 

WT^KG^Uue^WDiff = n^,,,,])! 2 < oo, (26) 

[l2,C2] 
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ll#ifG,jr; n n [JI;Cl] || D ( // = ^ \QKc,Kin(n[ S2 , C2] ,U [suCi] )\ 2 < oo. (27) 



Given a diffeomorphism invariant spin-scalar network function II[, ljCl ], there are only 
finite number of terms H[j liCl ][W^ c) n s>C (=[ JiC ]] which are nonzero in the sum over equiv- 
alent classes [s, c] in (I22> . On the other hand, given one spin-scalar network func- 
tion n,, ee [ J C ], there are also only finite number of possible II[ J2]C2 ] such that the terms 

n[s 2 ,c2][W^ c )n SjCe [> (C ]] are nonzero. As a result, only finite number of terms survive in 

both sums over [52, C2] in Eqs. d26i and ( 1271 . Hence both ^Hkg.ip an d ^KCKin are well 
defined. Then it follows from Eqs. ( 1221 . J24I and d25i that they are positive and sym- 
metric operators densely defined in fiDiff, whose quadratic forms coincide with Qkc,4> 
and Q.KG,Kin on their form domains. Hence both Q.kg$ an d QiccKm have positive clo- 
sures and uniquely induce self-adjoint (Friedrichs) extensions of 'Hkcq an d 'Hkgxux 
respectively |28|, which we denote by c Hkg4 an d "J^kckw as well. As a result, the 
Hamiltonian operator 'Hkg defined by Ea. (l23> is also positive and symmetric. Hence 
it has a unique self-adjoint (Friedrichs) extension. 
□ 

We notice that, from a different perspective, one can construct the same Hamilto- 
nian operator 'Hkg without introducing an inner product on T>* . The construction is 
sketched as follows. Using the two well-defined operators W^'f and a s in d!4i . as 

well as their adjoint operators (W^'f)' and (W^)', one may define two operators on 
'Hmff corresponding to the two terms in fl!2i by 

w KG , Kin v Diff )[f 7 ] = e j^w Diff[ ^ x s J M 6 (Ky)6 n^nUi 

(28) 

In analogy with the discussion in section 5 and Ref. 11181 . it can be shown that both 
above operators leave 'Hoiff invariant and are densely defined on 'Hoiff- Moreover, the 
quadratic forms associated with them coincide with the quadratic forms in <I22I . Thus 
the Hamiltonian operator 'Hkg '■ - 'Hkg,^ + ^RCXm coincides with the one constructed 
in the quadratic form approach. 

In summary, we have constructed a positive self-adjoint Hamiltonian operator on 
'Hoiff for the polymer-like scalar field, depending on a chosen parameter A. Thus there 
is an 1 -parameter ambiguity in the construction. However, there is no UV divergence in 
this quantum Hamiltonian without renormalization, since quantum gravity is presented 
as a natural regulator for the polymer-like scalar field. 

4 Quantum Hamiltonian constraint equation 

In this section we consider the whole dynamical system of scalar field coupled to grav- 
ity. Recall that in perturbative quantum field theory in curved spacetime, the definition 
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of some basic physical quantities, such as the expectation value of energy-momentum, 
is ambiguous and it is extremely difficult to calculate the back-reaction of quantum 
fields to the background spacetime |29|. This is reflected by the fact that the semi- 
classical Einstein equation, 

R#[g] ~ \R\g\gafi = K< tafilg] >, (29) 

is inconsistent and ambiguous |30||1|. One could speculate on that the difficulty is 
related to the fact that the present formulation of quantum field theories are back- 
ground dependent. According to this speculation, if the quantization programme is 
by construction non-perturbative and background independent, it is possible to solve 
the problems fundamentally. In loop quantum gravity, there is no assumption of a pri- 
ori background metric at all. The quantum geometry and quantum matter fields are 
coupled and fluctuating naturally with respect to each other on a common manifold. 
On the other hand, there exists the "time problem" in quantum theory of pure gravity, 
since all the physical states have to satisfy certain version of quantum Wheeler-DeWitt 
constraint equation. However, the situation would be improved when matter field is 
coupled to gravity. In the following construction, we impose the quantum Hamiltonian 
constraint on < Hkin, and thus define a quantum Wheeler-DeWitt constraint equation for 
the scalar field coupled to gravity. Then one can gain an insight into the problem of 
time from the coupled equation, and the back-reaction of the quantum scalar field is 
included in the framework of loop quantum gravity. 

Recall that the gravitational Hamiltonian constraint operator ^Hq^N) can be well 
defined in by the uniform Rovelli-Smolin topology 1 14||5 1. Hence it is also well 
defined in the coupled kinematical Hilbert space < Hkm- Its regulated version via a state- 
dependent triangulation T(e) reads 



<H e GH (N) = <H e E (N)- 2(1 +y 2 )t%N) 
16 



e ijk x 



J Y v eV(a) KA)=v 

Tr(A( ff , 7 (A)r '^(A))- 1 [A(s k (A)), V uf ]) P X 



E(v) 



Tr^toCA))" 1 [A( Si (A)\ K E ]A(sj(A)T' [A( Sj (A)), fc] x 

A(5 i (A))- I [A(^(A)),t> t/ .])^-. (30) 
E{v) 

We now define an operator in Hian corresponding to the scalar field part fixciN) of the 
total Hamiltonian constraint functional, which can be read out from Eqs. Q and @ as 

where 

?WA0 = f d i xN—^6<jp<!P b j (d a <p)d b <!>, 

2 J 2 ^J\de^[\ 
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In analogy with the regularization and quantization in the previous section, the regu- 
lated version of quantum Hamiltonian constraint ( H'j CG (N) of scalar field is expressed 
via a state-dependent triangulation T(e) as: 

i>eV(y) 

where the operators 



Pa te,A Pa 

v(A)=v 





v(A)=v V^OO V^OO v(A)= 

are all cylindrically consistent up to diffeomorphisms. Hence the family of Hamiltonian 
constraint operators OH is also cylindrically consistent up to diffeomorphisms, and the 
regulator e can be removed via the uniform Rovelli-Smollin topology, or equivalently 
the limit operator dually acts on diffeomorphism invariant states as 

(^ c MT w )[f r ] = limV/[^ G JWfyl (32) 

for any f y e Cyl y (^ C ){3l) ® Cyl y ( StC )( f lf), Similar to the dual of "HcRiN), the operator 
< H' KG (N) fails to commute with the dual of finite diffeomorphism transformation oper- 
ators, unless the smearing function N(x) is a constant function over S. Note that the 
diffeomorphism invariant Hamiltonian operator "Hkc defined in the previous section 
is actually 9-(' KC (l). From Eg. OH . it is not difficult to prove that for positive N(x) 
the Hamiltonian constraint operator < Hkg(N) of scalar field is positive and symmetric 
in Hun and hence has a unique self-adjoint extension. It is pointed out in Ref.|31 1 
that, there can be problems associated with symmetric constraint operators for systems 
where the constraints close with structure functions as is the present case. However, 
not all the assumptions underlying this conclusion are valid in the framework of loop 
quantum gravity. For example, it is assumed in Ref.|31 1 that all the classical canonical 
variables and constraints could be promoted as well-defined operators in the kinemat- 
ical Hilbert space. However, it is well known that the classical connection and diffeo- 
morphsim constraint cannot be represented as well-defined operators in loop quantum 
gravity. This issue related to the symmetric Hamiltonian constraint operator was fully 
discussed in Ref.| 14). 

Our construction of < Hkg(N) is similar to that of the Higgs field Hamiltonian con- 
straint in Ref . 1 22 1 . However, like the case of 'Hug, there is an 1-parameter ambiguity 
in our construction of 'H/ i:G (A0 due to the real scalar field, which is manifested as the 
continuous parameter A in the expression of in J15I . Note that now A is not re- 
quired to be a constant, i.e., its value can be changed from one point to another. This 
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issue of ambiguity will be discussed again in a future publication |27|. Thus the total 
Hamiltonian constraint operator of scalar field coupled to gravity has been obtained as 

<H(N) = <H G r{N) + <H KC {N). (33) 

Again, there is no UV divergence in this quantum Hamiltonian constraint. Recall that, 
in standard quantum field theory the UV divergence can only be cured by renormal- 
ization procedure, in which one has to multiply the Hamiltonian by a suitable power 
of the regulating parameter e artificially. While, now e has naturally disappeared from 
the expressions of J23l > and fl33l >. So renormalization is not needed for the polymer-like 
scalar field coupled to gravity, since quantum gravity has been presented as a natural 
regulator. Together with the result in the previous section, this heightens our confidence 
that the issue of divergence in quantum fields theory can be cured in the framework of 
loop quantum gravity. The desired matter-coupled Wheeler-DeWitt equation can be 
well imposed as: 

- (n' KG m¥m ff )Vy\ = {rt CR {N)v Diff )[f y ]. (34) 

Note that the scalar field part < Hkg(N) acts nontrivially on gravitational quantum states. 
This can be regarded as the "back-reaction" of quantum matter field to the quan- 
tum gravitational field. On the other hand, comparing Ea. J34> with the well-known 
Schordinger equation for a particle, 

d d 
ih—\p{x, i) — H(x, —ih-—)\p{x, t), 
ot ox 

where ifr(x, t) e L 2 (R, dx) and t is a parameter labelling time evolution, one may take 
the viewpoint that the matter field constraint operator < n' KC {N) plays the role of ifij r 
Then (f> appears as the parameter labelling the evolution of the gravitational field state. 
In the reverse viewpoint, gravitational field would become the parameter labelling the 
evolution of the quantum matter field. 

5 Master constraint programme 

In order to avoid possible quantum anomaly and find the physical Hilbert space of 
quantum gravity, master constraint programme was first introduced by Thiemann in 
1161 . The central idea is to construct an alternative classical constraint algebra, giving 
the same constraint phase space, which is a Lie algebra (no structure function) and 
where the subalgebra of diffeomorphism constraints forms an ideal. Self-adjoint master 
constraint operators for loop quantum gravity are then proposed in Refs. 1171 and 1 18 1. 
The master constraint programme can be generalized to matter fields coupled to gravity 
in a straightforward way. We now take the massless real scalar field to demonstrate 
the construction of a master constraint operator according to the strategy in Ref.| 18 1. 
By this approach one not only avoids possible quantum anomaly which might appear 
in the conventional canonical quantization method, but also might give a qualitative 
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description of the physical Hilbert space for the coupled system. We introduce the 
master constraint for the scalar field coupled to gravity as 



3 \C(x 
r — ;^^=^— 



M := - d\—-===, (35) 



where C(x) is the Hamiltonian constraint in (|4}. After solving the Gaussian constraint, 
one gets the master constraint algebra as a Lie algebra: 

{<V0), <V(rf')} = 'V([$, $']), 
{V(N), M) = 0, 

{M, M) = 0, (36) 

where the subalgebra of diffeomorphism constraints forms an ideal. So it is possible 
to define a corresponding master constraint operator on *%>;//■ In the following, the 
positivity and the diffeomorphism invariance of M will be working together properly 
and provide us with powerful functional analytic tools in the quantization procedure. 

The regulated version of the master constraint can be expressed via a point-splitting 
strategy as: 



M< := - | <l\ | d'xx e (x-y) C g_ C g_ . (37) 



Introducing a partition f of the 3-manifold S into cells C, we have an operator H[ 

acting on any cylindrical funct 
gulation T(e), 

PA feA Pa 



_ _ C,y 

y G Cyl y (^l) ® Cyly(1i) via a state-dependent trian- 



% = Z **v) z -JM 



E(v)" CR ' V ^fW) 



veV(y) v(A)=v 

Z ^c(v)[^(^;;,.) 1 w; : ; . + (w^,.„) 1 w^,,], (38) 

veV(y) 



where 



hf R , v = j^e ijk Tv[A( aij (A)r'A(s k (A)r l [A(s k (A)), ^V uf ]) 
+ (1 + y 2 ) 3 .^^ 3 ^Tr(A(.,(A))-' [A( i( (A)), £ e ] 



x A(^(A))- 1 [A( ij (A)),^]A(^(A))- 1 [A(^(A)), ^]), 

" , ( A)=v V^(v) 4 
x Tr(r iA(i m (A)) _1 [A(s m (A)), ^^(^(A))" 1 [A(s„(A)), V*{ 8 ]) A> x 



£A js.(„\Jmn 
v fct v y/W) 



- t.Y I /'A i 
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x Tr^KA))- 1 lA( Sl (A)), ^^(^(A))- 1 [A(s m (A)% V^ 12 ] 

x A( i „(A))- 1 [A( i „(A)),y^ e 12 ])-^. (39) 

The notations here are as same as those in section 3. Note that Hi is similar to the 

Hamiltonian constraint operator "Nil) defined in last section, but is now divided by 
the square root of volume operator. Hence the action of H e r on a cylindrical func- 
tion f y adds analytical arcs a,y(A) with j -representation (or arbitrary chosen spin j- 
representation of S U (2) if one uses non-fundamental representations to express the 
holonomies in (1301 ') and points at f(s,-(A)) with representation A with respect to each 
vertex v(A) of y. Thus, for each e > 0, Hi, is a S U(2) gauge invariant and diffeomor- 

phism covariant operator defined on Cyl y (jK) ® CylyiU). The family of such operators 
with respect to different graphs are cylindrically consistent up to diffeomorphisms. So 
the inductive limit operator He is densely defined on Ukui by the uniform Rovelli- 
Smolin topology. Moreover, the adjoint operators of H e r , which are also cylindrically 
consistent up to diffeomorphisms, read 

<%)' - I «w Z -j=<0' j= 

+ X XcivWW-;/^ . + (40) 

veV(y) 

The inductive limit operator of (I40i is denoted by (#£) ' , which is adjoint to He as 

< gy ' > #c/y >toi = <?yS #C,y/y >kin = < (Hc,y) ^gy',fy >kin 

= < (Hcfgy'Jy >km=< (ficf 7 ,gy,fy >km • (41) 

Then a master constraint operator, M, in 'Hoiff can be defined as: 

(MT Di// )[n s , c ] lim Y D ///[V l -H c (H c yil s , c ]. (42) 

CeP 

Since ffg(«^)^II, jC is a finite linear combination of spin-scalar-network functions on 
an extended graph with skeleton y, the value of (1VPFd;//)[IIj iC ] is finite for a given 
^Diff £ ^Diff- So M^oiff is in the algebraic dual of the space of cylindrical functions. 
Moreover, we can show that it is diffeomorphism invariant. For any diffeomorphism 
transformation <p, 

0'^ mff )[frl = Mm \^ c 0p%f 7 ] 

~* A "* CeP 

CeP 

lim V Diff [ V i^(^)Vy], (43) 



CeP 
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for any cylindrical function f y , where in the last step, we used the fact that the diffeo- 
morphism transformation <p leaves the partition invariant in the limit P — > E and relabel 
<p(C) to be C. So we have the result 

{(r$W Di f f )Uy\ = (MV Diff )[f 7 ]. (44) 

On the other hand, given any diffeomorphism invariant spin-scalar-network state Tl[ S ,c], 
the norm of the result state MII[j ]C ] can be expressed as: 

\\Mn ls j\ Diff 
= ^ i < tfm [JiC] |n [JV] > Diff i 2 

[.s'.c'] 

[s',c'] ~* ' £ ' e ~* CeP 

= y Hm _J_ y y 

[s',C] WW <peDiff/Diff rM <f>eGS yM 

x < f/^f/^n JiCe [ JjC ]| ^ -i7 c (H c )'n sV < e [. sV <] > K i„ | 2 

= y, Hm _J_ y y 

[*'■<"] rtI,t) ipeDiff/Diff^ <p'eGS yM 

s,ce[s,e]|ris',c'e[s',e'] >A'in I 

CeP 

= Z UJ^^o^'^C^c^c) 1 ^,^]]! 2 , (45) 

[s',c'] ' £ ' £ ^ CeP 

where we make use of the fact that M commutes with diffeomorphism transformations. 
Note that, the cylindrical function YiCeP \H^.(H^,yTl sce [ StC ] is a finite linear combina- 
tion of spin-scalar-network functions on some extended graph, so that there are only 
finite number of terms contributing to the sum in Ea.J45>. Hence it automatically con- 
verges. So the master constraint operator M defined by Eq. J42t is densely defined on 

^Diff- 

We now compute the matrix elements of M. Given two diffeomorphism invariant 
spin-scalar-network functions II[ JliCl ] and ri[ j2jC2 ], the matrix element of M is calculated 

as 

< n^^jiMin^c,] > Diff 
= (Mn [S2 , C2] )[n Sl , ClE[Jl , Cl] ] 

= P J™^Z l n[l!fl][ ^ )tn '" c,£[,1 " l] 

~* ,e ' e ^ CeP 

= i im y l —L- y y 

CeP y^-> <peDiff/Diff j( „, C2) <p'eGS r( „, C2) 
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x < JJ v U V 'Tl s ^ C2e [ s ^ C2 ]\H^{H e c )^Tl SuC2e [ SuC ^ > Kl 

= y um y 1 -— y y 



x < c/^f/^'nj 2jC2e [s 2iC2 ]|H^ r(s ^n J>c >Kin< n JjC |(^) 1 'iij liCl€ [j liCl ] >#,■„ 

= y y * 

[>,c] v€V(y(i,c€[j,c])) 

x n [S2>e2] [//^ (Sje) n !>ce[s , c] ] < n,^|(^)tn, ljCl6[JljCl] > ffin , (46) 

s,ce[.s,c] 

where Diff y is the set of diffeomorphisms leaving the colored graph y invariant, GS y 
denotes the graph symmetry quotient group Diff y /TDiff y where TDiff y is the dif- 
feomorphisms which is trivial on the graph y, and n y is the number of elements in 
GS y . Note that we have used the resolution of identity trick in the fourth step. Since 
only finite number of terms in the sum over spin-scalar-networks (s,c), cells C e 
f, and diffeomorphism transformations ip are non-zero respectively, we can inter- 
change the sums and the limit. In the fifth step, we take the limit C — > v and split 
the sum 2. s , e mto D[>,c] Zj,ce[>,c]> where [s, c] denotes the diffeomorphism equivalent 
class associated with (s, c). Here we also use the fact that, given y(s, c) and y(s',c') 
which are different up to a diffeomorphism transformation, there is always a diffeo- 
morphism if transforming the graph associated with H* cy ^s,c ( v e y( s > c )) to that of 

K>y( S >,c>) n s',c' ( v ' 6 rOV)) with <p(v) = V, hence U [s ^ C2] [H^ sc) U SrC€M ] is constant 
for different (s, c) e [s, c]. 

Since the term 2 s , ce[ . s , e ] < U s J(6^U Sl 

,ci€[si,a] >Kin is independent of the param- 
eter e', one can see that by fixing an arbitrary state-dependent triangulation T(e'), 



y <^s,c\0v)^si,cMsi,ci] >Kin 
s,ce[s,c] 

= 2j < U<ptls,c\(H$ )'n.S,,C 1 £[.!, ,Ci] >B» 

<p 

= t/ ^^-l( ( y), ) y (. s , e) n -'.c|n. Sl , ei 6[. Sl , Cl ] 



- n [sliC1 ][^ eV , ( (r(!c))r(j>c) n 5>c ], (47) 

where ^ are the diffeomorphism transformations spanning the diffeomorphism equiva- 
lent class [s, c]. Note that the kinematical inner product in above sum is non-vanishing 
if and only if tp(y(s, c))) coincides with the extended graph obtained from certain skele- 
ton y(si,ci) by the action of (Hf : )' and v e V(np(y(s,c))), i.e., the scale tp^ 1 (e') of the 
diffeomorphism images of the tetrahedrons added by the action coincides with the scale 
of certain tetrahedrons in y(s, c) and <p~ l (v) is a vertex in y(s, c). Then we can express 
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the matrix elements (I46> as: 

< n [il)Cl] |#i|ii te>C2] > Diff 

- Z Z 2 e ^ n ^.«][^v,r(i, c ) n *.' ; e[s.c]] n [si,ci][^ y ( J)C )n s , re [i, c ]] 

[s,c] veV(y{s,ce[s,c])) 

= Z Z 5(A'n [ . !2 , C2 ])[n. s , ce[ . ! , c] ](H,',n [J1 , fl] )[n J , ee[J , c] ]. (48) 

[s,c] vEV(y(i,ee[s,c])) 

From Eq.j48> and the result that the master constraint operator M is densely defined 
on 'Hoiff, it is obvious that M is a positive and symmetric operator on 'Hoiff- Hence, 

it is associated with a unique self-adjoint operator M, called the Friedrichs extension 
of M. We relabel M to be M for simplicity. In conclusion, there exists a positive and 
self-adjoint operator M on tiDiff corresponding to the master constraint ( I35> . It is 
then possible to obtain the physical Hilbert space of the coupled system by the direct 
integral decomposition of 'Hmff with respect to M. 

Note that the quantum constraint algebra can be easily checked to be anomaly free. 
Ea. J44t assures that the master constraint operator commutes with finite diffeomor- 
phism transformations, i.e., 

[M, f>;] = 0. (49) 

Also it is obvious that the master constraint operator commutes with itself, 

[M,M]=0. (50) 

So the quantum constraint algebra is precisely consistent with the classical constraint 
algebra ( I36> in this sense. As a result, the difficulty of the original Hamiltonian con- 
straint algebra can be avoided by introducing the master constraint algebra, due to the 
Lie algebra structure of the latter. 
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